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Abstract. 

The idealized mode coupling theory (MCT) is applied to colloidal systems 
interacting via short-range attractive interactions of Yukawa form. At low 
temperatures MGT predicts a slowing down of the local dynamics and ergodicity 
breaking transitions. The nonergodicity transitions share many features with the 
colloidal gel transition, and are proposed to be the source of gelation in colloidal 
systems. Previous calculations of the phase diagram are complemented with 
additional data for shorter ranges of the attractive interaction, showing that the 
path of the nonergodicity transition line is then unimpeded by the gas-liquid critical 
curve at low temperatures. Particular attention is given to the critical nonergodicity 
parameters, motivated by recent experimental measurements. An asymptotic model 
is developed, valid for dilute systems of spheres interacting via strong short-range 
attractions, and is shown to capture all aspects of the low temperature MCT 
nonergodicity transitions. 
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1. Introduction 



Colloidal particles interacting with strong attractions aggregate to form interesting 
structures. One such instance is gel formation, which occurs when the interparticle 
attraction is strong and of short range relative to the particle size. Gelation of 
colloidal systems appears to be a common, if not universal, phenomenon as it has been 
observed experimentally in numerous quite different colloidal-like systems. Notable 
examples of gel forming systems include colloid-polymer mixtures ^ ^, |^, |[ ^, |^ , in 
which the polymer is non-adsorbing and of low molecular weight, and suspensions of 
stcrically stabilized particles in marginal solvents j|, ^, O^FsL Emulsions , 
emulsion-polymer mixtures J^5|, colloid-surfactant mixtures |l6|, |l7|, and globular 
protein systems |l9[ 20, |2l[| are other examples of colloidal systems that exhibit 
gel formation. 

The phase diagrams of the colloid-polymer mixtures have been examined in detail, 
revealing that they are of the gas-solid type without a triple point and without a 
liquid phase when the attraction is short ranged. The disappearance of the liquid 
phase is a well-documented effect |2^, |^ |2^, |6[ |2^, |8|, which is caused by 
the restricted range of the attraction. This situation is rather unique to colloidal 
systems in that most molecular attractions are of comparable range to the molecular 
dimension, precluding such phase behavior. 

Whereas the equilibrium phase behavior of these systems is well understood, 
a fundamental understanding of the gel transition has been more difficult to 
achieve. Colloidal gels are characterized by ramified structures with particles located 
predominantly in clusters |6[ 0. It is thus natural to attempt to describe this 
phenomenon with percolation theories. This may seem especially appropriate in 
that density fluctuations are very slow near the gel transition, and particle clusters 
behave as nearly static objects [ pT| . However, rather poor agreement results when the 
percolation transition is compared to the experimental results for the gel transition 
line in the phase diagram, showing that the density dependence of the percolation 
transition is too strong |l^ . 

Another suggestion has been forwarded in which the gel transition is attributed 
to dynamic percolation within a gas-liquid binodal which is metastable with respect 
to gas-solid coexistence 0. Such metastable binodals have been observed also in 
solutions of globular proteins. However, gel (or precipitate) formation appears to 
occur in these systems also in regions of the phase diagram outside the metastable 
binodal |Q , suggesting that the gel transition is not only triggered by an instability 
towards local density fluctuations. The same seems to hold for sterically stabilized 
suspensions as well, which exhibit a gel transition at supercritical temperatures pj]] . 
Nevertheless, the same suspensions do form gels also under metastable conditions 

, as do colloid-polymer mixtures ^, ||, 0, |9| , prompting careful studies of any 
connection between the gel transition and metastability |0, |2^. 

We have recently suggested an explanation for the gel transition in Ref. jSOt , 
referred to henceforth as I. In this scenario colloidal systems form gels as a result of 
an arrested structural relaxation due to the self-trapping mechanism called the cage 



effect |31|, |3^, |3^, 34 , which is the same effect often thought to be responsible for 
the the liquid-glass transition. The idealized mode coupling theory (MCT), aimed at 
describing the cage effect in dense liquids, was seen to contain a bifurcation separating 
ergodic from nonergodic motion also for systems of particles interacting via strong 
short-range attractive interactions. In I we attributed the nonergodic states to gel 
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formation; hence, our suggestion is that the gel transition can be described within the 
same theoretical framework as the liquid-glass transition. 

Many of the qualitative aspects of the gel transition were found to be reproduced 
by the MCT calculations for the hard core attractive Yukawa system (HCAY). The 
calculated phase diagrams were found to exhibit gel transition lines that connect 
smoothly with the hard sphere glass transition at high temperature and extend to the 
critical and subcritical regions at low temperature along paths that depend critically 
on the range of the attraction. The phase diagram obtained by Verduin and Dhont 
is qualitatively reproduced by the MCT when the attraction is of intermediate range 
such that the gel transition meets the critical point. For shorter-range attractions 
the gel transition passes above the critical point, suggesting that structural arrest 
occurs instead of gas-liquid phase separation, which appears to agree with some 
measurements on sterically stablized particle systems |l^, |l^, |l^ . 

A finite zero-frequency shear modulus is predicted by MCT, in agreement with 
many measurements on colloidal gels 0, 10, 1^, which is expected to be intimately 
connected to a finite yield stress observed in other measurements At present 

the theory cannot account for the density dependence of the modulus , indicating 
that the mesoscopic structure of the gels is important and is not captured correctly. 
In addition, the theory alone cannot explain the growth of the small angle peak of 
the static structure factor nor the fractal scaling of its peak position observed on 
quenching of many suspensions It does, however, provide an explanation for 
why the growth of the small angle peak slows down and arrests following deeper 
quenches. Very differently from earlier approaches, which attributed the ultimate gel 
arrest to a global jamming or percolation transition, our microscopic theory predicts 
the structural arrest to be driven by anomalies of and a slowing down of the local 
dynamics. 

The purpose of this article is to provide a more detailed account of the results of 
the model study in I. In particular, the nonergodicity parameters (cf. Section 2.1) 
for several attraction ranges will be shown, motivated by the recent measurements 
on colloid-polymer systems by Poon and co-workers [Q. In addition, a detailed 
description is given of the asymptotic model developed in I. This model captures the 
relevant features associated with our suggested scenario for gel formation. A discussion 
of the relevance of these results is also included. 



2. Theory of colloidal gelation 

2.1. Mode coupling theory 

The idealized mode coupling theory (MCT) assumes that the dominant mechanism 
for structural relaxation in dense liquids is the cage effect. At short times particles 
are trapped in the surrounding cage of neighboring particles. At longer times 
particle escape from the cage leads to structural relaxation to equilibrium. For 
sufhciently strong interactions a bifurcation occurs in the governing equations — 
interpreted as the permanent entrapment of particles within their cages — causing the 
intermediate scattering function Fq (t) to acquire a non-zero long-time limit known as 
the nonergodicity parameter, glass form factor, or Debye- Waller factor fq ~ Fq{t ^ 
oo) /Sq, where Sq is the structure factor and q is the modulus of the wavevector. The 
transition from a vanishing to a finite long-time limit of Fq(t) is discontinuous and 
defines the liquid-glass transition within MCT. It is not a conventional thermodynamic 
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phase transition, but rather an entirely dynamic transition, interpreted generally as an 
ergodic-nonergodic transition. This simplified scenario captures many aspects of the 
liquid-glass transition in molecular liquids as well as hard sphere suspensions 

|2|,||,0||,|^,0. 

The governing MCT equation for the time evolution of the intermediate scattering 
function reduces in the long-time limit to 



■J dkF(q,k)25,5fc5|q_k|/fc/|q-k| 



1 - /, 2(2^)3g2 

^(q,k) =q- (q-k) C|q_j^l -q-kcfc (1) 

where p is the number density and Cq — (1 — S^^)/p, which appears in the vertex 
function T^(q, k), is the Fourier-transformed direct correlation function. Consideration 
of the single particle motion leads to another set of equations for the incoherent 
nonergodicity parameter (also known as the Lamb-Mossbauer factor) 



fs 
Jq 



j dkF^(q,k)25fc/fe/|^q_]^| 



1 - f^q (27r)3g2 
F^(q,k) -q-kcfc (2) 

where — F^{t ~^ oo), with the self-intermediate scattering function F^{t). 

For a given Sq, equations 1 and 2 are closed equations for fq and /|. They are 
solved numerically by iteration starting from the initial iterate fq — fq — 1. With 
this starting point the iteration converges monotonically to the correct solution for 
the nonergodicity parameters, given by the largest solution to the equations The 
wavevector integrations are performed efficiently using Simpson's rule on a uniformly 
discretized grid: q — iAq, i — 0, . . . , N. The critical glass transition boundary is 
identified by bracketing of the given input conditions, such as temperature and density, 
which delineates regions where fq is zero and finite; regions of the phase diagram 
which result in finite nonergodicity parameters are identified as glass states. Most 
calculations were done with the parameters Aq — 0.3(t~^ and N = 800. 



2.2. Model system 

The hard core attractive Yukawa (HCAY) interaction potential captures both short- 
range excluded volume interactions and variable-range attractive particle interactions. 
It is given by 

u{r)/kBT ^ I l,(^r/a-l) !!<!.'^'^ (3) 

L r/cr 

where the dimensionless parameter K regulates the depth of the attractive well and 
the reduced screening parameter b sets the range of the attraction. 

The equilibrium phase diagrams for a number of attraction ranges have been 
obtained by computer simulations p7| , showing indeed that the liquid phase 
disappears upon restricting the range of the attraction. Because the HCAY static 
structure factor is known semi-analytically from the mean spherical approximation 
(MSA) ^ Q, it is convenient to base our study on this case. The penalty is 
that the resulting predictions are limited not only by the approximations made in the 
MCT, but also those made in the MSA. Therefore, the results should be viewed as 
qualitative rather than quantitative. 
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2.3. Asymptotic model 

In the following we describe an asymptotic model, originally developed in I, that 
captures the relevant features of the suggested gelation mechanism. At low densities 
((/) = ■Kpa^/Q 0) the Ornstein-Zernike direct correlation function becomes 
independent of density. Specifying this to the MSA of the HCAY fluid, we obtain 
Kb{a/bfd{q,a,b), where c(g,cr,&) = (b/a)"^ J^^^dre-''^''/'^-^'>e''^-^/r in which 



we have neglected the contribution from the hard core which is independent of K and 
b at low densities. In the limit of strong attractive interactions at low densities, the 
following scaling simplifies the MCT equations 

(j) ^ Q and oo, so that F = — ^ — = constant (4) 

In addition, 5*^ — )■ 1 in this limit and the MCT equations (1 and 2) simplify 
considerably such that the nonergodicity transitions occur at F = Fc(fe), leading to 
the asymptotic prediction cx 1/ -y/^- 

For short-range attractive interactions, in the limit & — > oo, a further simplification 
arises because the MCT vertex functions become linear functions of F only, if we 
introduce the rescaled wavevectors q ~ qa jb. Equation 1 simplifies to 



^ = -ydky(q,k)V,/„_fc, (5) 

with the dominant contribution to the rescaled vertex function as 
l/(q,k)2 = V^^(q,q-k)2 + F«(q,k)2 

y-(q,k)^^A jq-^)', (6) 

where the nonergodicity parameters depend only on the rescaled wavevectors fq-^fq- 
The same set of equations is found to govern the incoherent nonergodicity parameters 
/|. Note that in this limit single particle and collective nonergodicity factors are 

identical and exhibit the small wavevector expansion fq = l — q^ (rgb/a)'^, where is 
the localization length, or root-mean-square displacement, in the glass. 

In this limit the coupling constant F assumes a unique value at the transitions, 
which is found from solving equations |^ and 6 numerically resulting in Fc « 3.02 for 
b ^ oo. We emphasize that these asymptotic forms are valid for the HCAY systems 
using the MSA in the prescribed limits. 



3. Results 



3.1. Phase diagrams 

Using the MSA static structure factor |4|, H, ||] as input, the MCT was solved for 
several screening parameters: b — 7.5, 20, 30, and 40. The progression of the glass 
transition can be traced from the (Percus-Yevick) hard sphere limit, corresponding 
to K=0, to lower temperatures in terms of the reduced temperature K~^. In I we 
showed the phase diagrams for b ~ 7.5, 20, and 30. This work adds an additional 
phase diagram for 6 = 40. 

As shown in I, at low temperatures (large values of if) the glass transition is traced 
along different paths in the phase diagrams depending on the value of the screening 
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parameter, i.e. the range of the attractive interaction. In all cases they bend towards 
lower densities when the temperature is decreased sufficiently. The transition lines 
for intermediate-range attractions (6 = 7.5 and 20) reach subcritical temperatures 
on the liquid side of the spinodal. For shorter ranges of the attraction {b — 30) the 
nonergodicity transition line lies entirely within the fluid phase above the two phase 
region, and extends to subcritical temperatures at low densities on the vapor side of 
the spinodal. This is shown in detail in figure in which the results for 6 = 30 and 
40 are shown. As seen, decreasing the attraction range further to 6 = 40 causes the 
nonergodicity transition line to move away from the spinodal curve. In contrast to 
the MSA phase diagrams studied in I, the 6 = 40 transition line is located sufficiently 
far away from the spinodal curve to remove the additional nonergodicity transition 
line that appears in the MSA phase diagrams with lower b. This type of transition is 
discussed in the Appendix of I. Also shown in figure are the spinodal curves, defined 
by the condition Sq oo for q — > 0. In the present context they are shown as an 
indication of where gas-liquid phase separation is likely to occur, provided a liquid 
phase is present. 

As pointed out in I, the MCT nonergodicity transition lines resemble the gel 
transition lines determined experimentally. Most notably, the diagram with b = 20 
(see I), which displays a near meeting of the nonergodicity transition and the critical 
point, bears a strong resemblance to the phase diagram determined by Verduin and 
Dhont [ pT| in their measurements of sterically stabilized silica suspensions. Other 
measurements on sterically stabilized suspensions found a gel transition line with no 
apparent evidence of phase separation, which can be expected, according to the MCT, 
for attractions of very short range, such as 6 = 30 or 40 in figure 0. 

One important aspect of the MCT nonergodicity transitions is that they are not 
induced by long-range structural correlations associated with critical fluctuations. 
This is demonstrated by the asymptotic model which is valid in the limit in equation 
^ and 6 oo for which Sq ~ I. The predictions of the model are shown in 
figure as chain curves. There is quantitative agreement at low densities and the 
model predictions reproduce the results of the full calculation for the nonergodicity 
transitions qualitatively at moderately high densities. The agreement with the 
asymptotic model demonstrates that the q limit of Sq is decoupled from the 
structural arrest at low densities and temperatures. This conclusion is of particular 
importance as not all relevant mode couplings expected for systems near the critical 
point [E3 are included in this version of MCT. 



3.2. Nonergodicity parameters 

Dynamic light scattering measurements of the nonergodicity parameters led, inter 
alia, to the identification of the glass transition for the colloidal hard sphere system 
p7| , m, |o| . As MCT predicts very different behavior for these when strong short- 
range attractions are present [pof , it can be expected that such measurements will 
provide a similarly decisive test of the MCT of gelation. For that reason we focus 
here on the behavior of the nonergodicity parameters as functions of density and 
temperature. 

When the attraction is of moderately short range, such as for b — 7.5, figure ^ 
shows that the coherent nonergodicity parameters deviate little from the hard sphere 
fq as the temperature is lowered, with the exception of an increasing q value. 
Note that the fq shown in figure |2| and later figures are the critical nonergodicity 
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parameters, which appear along the MCT transition hnes shown in figure 2 of I and 
figure of this work. The increase of the q value results in this case mainly from 
the increased contribution of long-range correlations, which derive from the increase 
of the q ^ limit of Sq upon approaching the liquid side of the spinodal curve. Shown 
also in figure |^ is the prediction from the asymptotic model based on equation ^ and 
& ^ oo, which does not capture the behavior of fq at the lowest temperatures. This 
should be expected as the attraction is not particularly short ranged. 

Decreasing the range of the attractive interaction to 6 = 20, which causes the 
nonergodicity transition line to move closer to the critical point (see figure 2 of I), leads 
now to form factors that deviate more from hard sphere behavior at low temperatures. 
Compared to 6 = 7.5, this shorter range attraction causes coupling among more 
wavevector-modes, leading to an increase in the width of fq as the temperature is 
decreased. This effect is produced by the increased width of the static structure 
factor, caused by strong short-range correlations due to the attraction. As in the case 
of 6 = 7.5, the g — s- limit of fq for b = 20 increases when the temperature is decreased. 
Now, however, the asymptotic model yields a reasonably accurate prediction for fq 
at the lowest temperatures, which can be expected from the rather close agreement 
between the asymptotic model and the full calculation of the transition line for 6 = 20. 

When the range of the attraction is decreased further by changing the HCAY 
screening parameter to & = 30, the path of the nonergodicity transition is unimpeded 
by the gas-liquid critical curve (see figure |l|) ; it extends to subcritical temperatures on 
the low density, vapor side of the spinodal curve. The form factors along the transition, 
shown in figure ^, are qualitatively similar to those corresponding to & = 20, except 
that the width of fq is further increased due to stronger short-range correlations. The 
agreement with the asymptotic model is improved at the lower temperatures where 
the oscillations in fq are now somewhat suppressed. The improved agreement should 
be expected because the asymptotic model is based in part on the limit — s- 0, which 
can be fulfilled by this system as the nonergodicty transition line lies entirely in the 
(very likely metastable) single phase fluid regime. 

The dynamics of the gel transitions can be expected to be anomalously stretched 
over many orders in time and to exhibit a rapid slowing down upon approaching the 
nonergodicity transition lines. As discussed in detail in reviews of the MCT |33| , 
this and a number of other qualitative aspects can be understood from the factorization 
property and asymptotic expansions, which describe the sensitive variation of the cage 
dynamics close to the transition. One finds: Fq{t)/Sq = fq + hqG^{t). As the so- 
called /3-correlator, G^{t), also exhibits numerous universal features depending on 
one material parameter only, the exponent parameter A, this expression provides for 
crucial experimental tests of MCT transition scenarios, such as those performed for 



hard sphere colloids |3§, |39[ We find A = 0.89 predicting a very anomalous 
stretching and rapid increase of the longest relaxation time |3^, ^ . Figure || shows 
the two wavevector dependent amplitudes, which describe the gel structure (/^, also 
included in figures ||- ^ and localized cage dynamics {hq). The amplitude factor hq, 
which describes the spatial extent of that dynamical process which arrests at the gel 
transition, is found to be peaked at rather large wavevectors stressing that the local 
motion of the colloids is suppressed at the gel transitions. 
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4. Discussion and Conclusions 

Structural arrest as described by the idealized MCT is found to be a plausible 
explanation for colloidal gelation. The underlying cause of gelation in this scenario is 
a breaking of ergodicity caused by strong short-range attractions in dilute systems. 
An ergodic-nonergodic transition is characteristic also for the glass transition within 
the framework of MCT. The low temperature gel transitions are accompanied by the 
cessation of hydrodynamic diffusion and the appearance of rather large finite elastic 
moduli due to the particles being tightly localized in ramified clusters (for more details 
see I). 

At relatively high temperatures, corresponding to weak attractions among 
the suspended particles, the cage surrounding a typical particle is distorted by 
dimerization (see also |Q ) . The cage must be reinforced by increasing the critical 
colloid density before structural arrest ensues, which leads to a glass transition line 
that moves initially towards higher density with decreasing temperature. This trend 
appears to require hard core repulsions and becomes more pronounced when the range 
of the attractive interaction decreases. It has been observed in the adhesive hard sphere 
system ^ and the Yukawa systems investigated in I and this study. For the 
colloid-polymer mixtures recrystallization of glassy samples has been reported upon 
introducing a small concentration of low molecular weight non-adsorbing polymer 
Q. We attribute this effect to the glass transition line possessing an initial slope in 
accordance with the MCT predictions, i.e. in the direction of higher density with 
decreasing temperature. 

Stronger short-range attractions cause particle aggregation, leading in effect to 
an increase of the density in the local environment of a typical particle. According 
to MCT, aggregation can lead to structural arrest of the long-time dynamics despite 
the bulk density being much lower than the hard sphere glass transition density. As 
shown in I and figure |l|, the precise path of the nonergodicity transition line in the 
phase diagram depends now critically on the range of the attraction, reminiscent of the 
dependence of the freezing line on the range of the attraction |2^, |2^, ^ ^ . 

The low density nonergodicity transitions, which appear only for sufhciently short 
ranges of the attraction, are not only caused by the excluded volume effect which 
dominates at higher densities, but are additionally affected by the low- and high-g 
behavior of the static structure factor. The asymptotic model demonstrates, however, 
that the short-range (high q) correlations become increasingly important as the range 
of attraction is decreased, such that they are the dominant cause of the structural 
arrest at low densities. This observation indicates that the nonergodicity transitions 
are unaffected by long-range correlations and the presence of the critical point. Thus 
the fact that the MCT transitions occur in close proximity to the spinodal curve for 
a significant range of b is merely coincidental. The results shown in figure |l| are in 
accord with this conclusion; the nonergodicity transition line for 6 = 40 is seen to be 
located further away from the critical curve than for & = 30. 

Nevertheless, several studies, particularly the extensive ones of the colloid-polymer 
mixtures, reveal what appears to be an intimate connection between the gel transition 
and the metastable gas-liquid binodal. It is possible that the gel transition changes 
character once it crosses the metastable binodal into the metastable and unstable 
regions. This is not in disagreement with our description of the arrest of the local 
dynamics, but indicates that the large-distance dynamics can be more complicated. 
Unfortunately we are yet unable to make any predictions inside the unstable region 
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with the present MCT because an appropriate Sq is not available and the theory 
assumes closeness to equilibrium [ p3| . 

Krall and Weitz have studied low density gels experimentally at small 
wavevectors, though still larger than the wavevectors characterizing the small angle 
scattering peak. The low density gels are associated with finite nonergodicity 
parameters, supporting the MCT gelation mechanism. However, Krall and Weitz 
measure relatively small values for the nonergodicity parameters, whereas the full 
MCT solutions, as well as the asymptotic model described in Section 2.3, predict large 
values for fq at small wavevectors. Such large values are observed experimentally also, 
but only at somewhat higher densities. We are attempting presently to locate a low 
density region in the MSA phase diagrams with this type of dynamics. 

Experimental measurements of the dynamics of gels at higher densities have been 
conducted. Sterically stabilized suspensions and colloid-polymer mixtures with 
low molecular weight polymer show that these gels are associated also with 
finite nonergodicity parameters, further supporting the present model calculations. 
Moreover, the measured fq assume values not too much below unity both at low q 
and at values near the q corresponding to the principal peak of Sq . 

Poon and co-workers ||2^ have made detailed studies of the low density and 
temperature region of the colloid-polymer phase diagram for low molecular weight 
polymers. They identified a variety of different dynamics which can be reconciled 
with the MCT calculations provided their system belongs to the HCAY diagram with 
b « 20. For such a situation the gel transition meets the spinodal on the liquid side, 
without interfering with the critical behavior along the vapor side. Conjecturing that 
the gel transition can exist in the unstable region as found experimentally by Verduin 
and Dhont [ pT| , we obtain a phase diagram qualitatively similar to that determined 
by Poon et al. We expect that quenching of a suspension at low to moderately 

low density results initially in normal phase separation dynamics, characterized by 
either spinodal decomposition or nucleation and growth. This process proceeds until 
the denser domains reach the critical density for gelation, subsequently arresting the 
structure, which would be consistent also with the measurements in Ref. [ pl| . Thus 
we anticipate that the so-called transient gelation region, discovered by Poon et al. 
p9| for colloid-polymer systems with short polymers, corresponds to a nonergodicity 
transition in the unstable region of the p hase diagram. 

The asymptotic model of Section 2.3 provides additional support for this 
interpretation. First, arrest of long-range structures is caused by an arrest of the 
local dynamics. In contrast to percolation approaches for the gel arrest, the local 
density fluctuations also exhibit slowing down and dynamical anomalies. Second, the 
tight localization of the particles provides a rationale for using concepts like bond- 
formation and sticking probabilities, even though the particles obey diffusive equations 
of motion. Third, the approach to unity of the collective nonergodicity or Debye- 
Waller factors for small wavevectors, f{q) 1 for g — > 0, indicates that the particles 
are bound to infinite clusters in this limit. Momentum conservation, or Newton's law 
of action-and-reaction, otherwise would require f{q ^ 0) < 1, as found for glasses. 

In summary, gelation in colloidal systems is attributed to a breaking of ergodicity, 
captured qualitatively by the MCT applied to the HCAY system. The calculated phase 
diagrams show many similarities with those determined experimentally for colloid- 
polymer mixtures and sterically stabilized suspensions. Decreasing the range of the 
attractive interaction sufficiently causes the gel transition to move further into the 
single phase region. The gross features of the MCT nonergodicity parameters at low 
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temperatures are seen to be in accord with the few existing measurements, supporting 
the proposed gelation mechanism. 
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Figure captions 



Figure 1. Enhancement of the low density and temperature region of the HCAY 
diagrams for b = 30 and 40. The MSA spinodal curves are shown with the critical 
points denoted by (•), together with the corresponding MCT nonergodicity transition 
lines as labeled. The chain curves correspond to the asymptotic prediction in equation 
^ and b —> oo with Fcib — > oo) = 3.02. The inset shows the phase diagram for b = 40 
for the same density and temperature ranges as in figure 2 of I. 



Figure 2. HCAY critical coherent nonergodicity parameters fq for 6=7.5 along the 
critical boundary as functions of normalized wavevector qa and Yukawa prefactor K. 
Note that the volume fraction varies according to the b=7.5 nonergodicity transition 
line in figure 2 of I. The K values, incremented by 1, run from bottom-to— top from 
to 7. The curve shown in bold is the asymptotic prediction resulting from equation 
fi and b — + oo . 



Figure 3. HCAY critical coherent nonergodicity parameters fq for 6 = 20 as in 
figure ^. The K values, incremented by 2, run from bottom— to-top from to 12. 
The asymptotic prediction from equation U is shown as the bold curve. 



Figure^ 4. HCAY critical coherent nonergodicity parameters fq for fe = 30 as in 
figure H. The K values, incremented by 2, run from bottom— to-top from to 16. 
The asymptotic prediction from equation U is shown as the bold curve. 
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Figure 5. Critical coherent nonergodicity parameter fq and amplitude factor hq as 
functions of the scaled wavevector for the asymptotic model in equations ^ - 6 and 
b oo. 



